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Abstract 

We investigate how entanglement spreads in time-dependent states of a 1-1-1 
dimensional conformal field theory (CFT). The results depend qualitatively on 
the value of the central charge. In rational CFTs, which have central charge below 
a critical value, entanglement entropy behaves as if correlations were carried by 
free quasiparticles. This leads to long-term memory effects, such as spikes in the 
mutual information of widely separated regions at late times. When the central 
charge is above the critical value, the quasiparticle picture fails. Assuming no 
extended symmetry algebra, any theory with c > 1 has diminished memory 
effects compared to the rational models. In holographic CFTs, with c ;§> 1, these 
memory effects are eliminated altogether at strong coupling, but reappear after 
the scrambling time t > /3 log c at weak coupling. 
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1 Introduction 

How does quantum information spread in a strongly coupled system, far from equi¬ 
librium? This question is central to the dynamics of a wide variety of experimental 
and theoretical systems, from cold atoms and condensed matter to quantum chaos and 
black holes (see for example [1], [2], [3] and [4,5]). 

Quantum information is encoded in entanglement, so a useful way to characterize 
its spread is through the time-dependent entanglement entropy of spatial subsystems. 
In this paper, we investigate this quantity in l-|-ld critical systems, in time-dependent 
states which have only short-range correlations at time t = 0. These states model a 
global quench from a gapped Hamiltonian. 

The entanglement entropy of a single, connected region in this context was com¬ 
puted by Galabrese and Gardy using conformal field theory (GET) [6-8]. In the limit 
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where the time t and size L of the interval are much larger than the initial correlation 
length the result is hxed universally by conformal symmetry, and depends only on 
the central charge c of the CFT. It is equal to c times the answer for a free boson. 
This can be modeled by assuming that entanglement is carried by free quasiparticles, 
propagating ballistically at the speed of light. 

The entanglement entropy of disjoint regions is a more detailed probe of how in¬ 
formation spreads. We focus on the case of two disjoint intervals A and B. In this 
case the mutual information I{A,B) = Sa + Sb — Saub bounds the strength of con¬ 
nected correlation functions of smeared local operators in the two regions [9]. Our 
main result is that the entanglement entropy SAusit) is not universal in general, but 
depends qualitatively on the value of the central charge: It interpolates between the free 
quasiparticle answer for small c, and the holographic prediction of three-dimensional 
quantum gravity at large c. 

Free quasiparticles vs. quantum gravity 

Suppose the regions A and B both have length L and are separated by a distance 
D > L, as in hgure la. In a theory of free quasiparticles, the entanglement entropy 
SAusit) is shown in hgure lb. After the quench, it grows linearly as entangled pairs 
spread, and saturates at the thermal value for temperature /3 ~ At t = D/2, 
entangled pairs created in the middle enter the two regions, as shown in hgure la, 
leading to a dip in Savjb and a corresponding spike in the mutual information. 

At the midpoint of the dip, the left-moving quasiparticles in region A are maximally 
entangled with the right-moving quasiparticles in region B, so Saub is exactly half the 
thermal value (after subtracting the divergent contribution at t = 0). In this scenario, 
maximally entangled degrees of freedom remain maximally entangled even as they 
propagate to very large separation. In other words, entanglement does not scramble. 

It was claimed in [6, 8] that the free quasiparticle behavior for multi-interval en¬ 
tanglement after a global quench is universal to all conformal held theories. This is 
surprising in a strongly coupled system, and we will argue that an assumption un¬ 
derlying this claim (about the analytic continuation of certain Euclidean correlators 
to Lorentzian signature) is justihed only in a restricted class of theories with a large 
degree of symmetry. A hrst hint of this was provided by holographic calculations of 
entanglement entropy [10-12], where, even in 1-1-1 dimensions, entanglement scrambles 
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Figure 1: (a) Quasiparticle picture for the entanglement entropy of two disjoint intervals. 
Entangled pairs are created at t = 0, which propagate in opposite directions at the speed of 
light. When one of the particles enters A and the other enters B, the entanglement entropy 
Saub decreases. (6) Memory effect in the entanglement entropy of two separated intervals 
after a global quantum quench. For intervals of size L separated by a distance D > L, the 
dip is centered at t = (D + L)/2 and extends over the range D/2 < t < D/2 + L. The solid 
line is the free quasiparticle answer, and the dashed line is the holographic answer. 


maximally — late time features in the entanglement entropy of widely separated sub¬ 
systems, such as the dip in hgure 16, are entirely absent. This apparent discrepancy 
was hrst observed in [10] and has also been explored in [13]. Our aim is to reconcile 
these two pictures, and to understand the middle ground. 

Results 

We argue that, in fact, the scrambling of entanglement is generic in l-|-ld conformal 
held theory: scrambling occurs whenever the central charge is greater than a critical 
value depending on the number of conserved currents. In a CFT without extended 
chiral symmetry, the stress tensor and its descendants are the only conserved currents, 
and entanglement scrambles if and only if c > 1. More generally, we dehne an effective 
central charge Ccurrents for the chiral sector of the theory, and hnd that entanglement 
scrambles when c > Ccurrents- 

The quasiparticle picture is accurate in theories with c = Ccurrents- These include 
the rational CFTs, as well as some non-rational theories like the compact boson at 
irrational radius squared. We refer to this class of theories as current dominated, for 
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the following reason. Denoting the density of states by Q{h,h), where {h,h) are the 
left and right conformal weights, c = Ccurrents implies 

D(h, h) D(h, 0) X D(0, h) as h,h^oo. (1-1) 

That is, the spectrum of high-dimension operators is dominated by conserved currents, 
which have dimension (h, 0) or (0, h). Correlators in such theories also factorize into 
left and right moving contributions, in a stronger sense than in general CFTs. 

The opposite extreme is a holographic CFT, which has Ccurrents ~ 1 and c 3> 1, 
and scrambles maximally. We recover this universal answer directly from our CFT 
analysis in the holographic limit. Theories with c > Ccurrents, but where c is not large, 
also scramble, but the entanglement entropy in these theories is not universal. These 
intermediate theories are not rational, so they cannot be solved exactly, and perhaps 
for this reason have received much less attention in the literature than rational CFTs. 
They are, however, perfectly normal, unitary quantum held theories, with a moderate 
central charge, a unique normalizable vacuum state, and a discrete spectrum when 
placed in hnite volume. Such CFTs with c > Ccurrents are much like their higher¬ 
dimensional cousins, since rational CFTs do not exist in d > 2. They have an inhnite 
but discrete set of primary helds. An example is a sigma model with a compact target 
space that has no isometries (or few isometries). Such a theory does not have the large 
quasiparticle dip shown in £g. 1, even in the scaling limit where times and lengths are 
much greater than the initial correlation length. An explicit example with c = 12 is 
discussed in section 5. 

The evidence for these claims is strong but not entirely conclusive. It follows from 
calculations of the Renyi entropies Sn in conformal held theory, which are related to the 
entanglement entropy by extrapolating n —)■ 1. We consider two models for a quantum 
quench: the boundary state model of Calabrese and Cardy [6-8] , and the thermal dou¬ 
ble model of Hartman and Maldacena [5]. In the thermal double model, we show that 
the S 2 Renyi entropy behaves precisely as claimed, and argue (somewhat incompletely) 
that the same applies for all Sn- Thus the evidence for scrambling in this model is very 
strong. In the boundary state model, our conclusions hold as long as boundary states 
in non-rational conformal held theory do not, for some unknown reason, produce a 
spurious singularity in correlation functions, beyond the usual singularities required 
by the operator product expansion, and the additional singularities from image points 
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across the boundary. This assumption is unproven but seems likely without evidence 
to the contrary, so we suspect that our conclusions hold also for boundary states. We 
therefore disagree with the claim in [6,8,14] that the quasiparticle picture for multiple- 
interval entanglement entropy after a quench is entirely universal — the derivation 
implicitly assumes that the theory is current dominated. The differences appear only 
for multiple intervals after a quench, or certain conhgurations of highly boosted inter¬ 
vals in vacuum, so the results for a single interval are unaffected. Similar differences 
would appear in the calculation of n-point correlation functions of local operators after 
a quench [15], but only for n > 2. 

We emphasize that our conclusions hold in the scaling limit where length and time 
scales are much larger than the correlation length in the initial state (as in [6-8]). 
The scaling limit is taken with c held hxed, and we assume c is large only when stated 
explicitly in section 4. In holographic theories, there is an additional consideration since 
we must take both the large-c limit and the scaling limit. This order of limits is not 
the origin of the discrepancy between holographic and rational CFTs. The interplay 
between the scaling limit and the large-c limit is addressed in section 4.4. At strong 
coupling, assuming a sparse spectrum of low-dimension operators, the conformal held 
theory answer reproduces the holographic prediction for mutual information after a 
global quench. In weakly coupled CFT, it leads to a quantitative prediction for stringy 
corrections to the holographic entanglement result. 

Outline 

The Calabrese-Cardy boundary state model for a global quench, and the related tech¬ 
nology of twist correlation functions in CFT, is reviewed in sections 2.1-2.2. In sec¬ 
tion 2.3, we show that around the time of the dip in the quasiparticle entanglement 
entropy, the behavior of the Renyi entropies is controlled by a light cone singularity 
in the twist correlator, and that this singularity is not hxed by the ordinary operator 
product expansion. 

The aim of section 3 is to understand the nature of this light cone singularity. 
First, we relate it to a technically simpler question about the mutual information of 
ohset intervals in a thermal double CFT. In this setting, in section 3.3 we compute the 
singularity in the S 2 Renyi entropy by conformally mapping the replica geometry to 
the torus partition function and give a dehnition of Ccurrents- Based on properties of the 
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spectrum in rational (or nearly rational) vs non-rational CFT, we conclude that the 
quasiparticle picture produces the correct S 2 only for central charge c below a critical 
value. In section 3.4 we attempt to extend this to general Renyi index n, but only 
partially succeed. Section 3.5 is devoted to the marginal case of the c = 1 free boson, 
and the question of whether the dip gets larger or smaller is addressed in section 3.6. 

In section 4 we explore the large c holographic limit. The gravity calculations of 
mutual information and second Renyi entropy after a global quench are sketched in 
sections 4.1-4.2; these are reproduced from CFT (with some extra assumptions) in 4.3; 
and the weakly coupled limit of a particular holographic CFT, the D1-D5 symmetric 
orbifold, is considered in section 4.4. 

In the conclusions, we briefly comment on several open questions and possible 
extensions of these results, including local quenches and entanglement negativity. 

2 Origin of entanglement memory 

Consider a H-ld CFT, in an initial pure state |4/) with a hnite correlation length ^ and 
hnite energy density. This state may be produced, for example, by globally quenching 
the Hamiltonian of a gapped system to a critical point. At time t = 0, regions of size 
T S> ^ exhibit area-law entanglement. The state is time dependent, and since it has 
hnite energy density, we may expect it to exhibit volume-law entanglement for a single 
interval at late times. This was beautifully demonstrated in [6-8]. 

A more detailed probe of entanglement in this system is the entanglement entropy 
or mutual information of separated intervals. Dehne the regions 

A : x G [0, L], B : X & [D + L, D + 2L] . (2.1) 

To simplify the discussion we will always assume 

t>0, D>L, (2.2) 

and only comment on the case D < L in the conclusions. 

The entanglement entropy Saub does not measure the entanglement of A with B, 
but rather the entanglement of A U R with the rest of the system — if A and B are 
highly entangled with each other, then they cannot be highly entangled with degrees 
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of freedom elsewhere, so Saub is small. It is UV divergent, but the divergence is 
time-independent, and could be removed entirely by considering instead the mutual 
information 

I{AB) = Sa + Sb-Saub . (2.3) 

Two possible behaviors of the entanglement entropy as a function of time are illustrated 
in figure 1. The linear rise at early times comes from correlations spreading until 
t ~ L/2. The dip, when it exists, can be understood as coming from entangled degrees 
of freedom that started at the midpoint between the two intervals, with left-movers 
entering region A and right-movers entering region B around t ^ D/2. When there is 
a dip in Saub, there is a corresponding peak in the mutual information. 

In this section we will setup the calculation of the entanglement entropy, following 
[6-8], and isolate the part of the computation responsible for the dip around t ~ D/2. 

2.1 Review of the boundary state model 

The experiment just described is in a particular pure state, |T). We cannot hope to 
capture the exact details of an arbitrary pure state, but it was argued by Calabrese 
and Cardy that the behavior of the entanglement entropy (and correlation functions) 
for f, L, D S> ^ is universal. The universal features of |T) are captured by the state 

1$) = (2.4) 

where \B) is a conformal boundary state, H is the CFT Hamiltonian, and /3 ~ The 
boundary state \B) is not normalizable, but evolution in Euclidean time spreads out 
correlations and renders the state normalizable. 

Entanglement entropy in the state |$) can be computed by the replica method. 
Define the Renyi entropy 

Sn = log Tr p" , (2.5) 

1 — n 

where p is the reduced density matrix of a subregion, and n > 2 is an integer. The 
entanglement entropy S = — trplogp is obtained by analytic continuation in n —)• 1. 
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Twist correlators 


In the state (2.4), the Renyi entropy can be computed by a path integral. For the 
region AU B dehned above, matrix elements of Pavjb at t = 0 are computed by the 
path integral on an inhnite strip: 


B 






{^'\Paub\p^) = \ 



\ ^ 


B 


( 2 . 6 ) 


where (p and p' denote states on region AVJ B. Reading this diagram from bottom to 
top, the boundary condition at the bottom dehnes the state |i3), the Euclidean time 
evolution over the bottom half of the strip produces and the path integral 

over the upper half corresponds to 

The replica partition function Tr p\jB is computed by the path integral on n copies 
of this system, glued along region AU B-. 


B B B 
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n' 

2 
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S' 
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B B B 

(2.7) 

with slits identihed as indicated by the numbers 1... n, V .. .n'. This is a path integral 
on a Riemann surface with n — 1 handles, and n holes where we impose boundary 
condition B. It can also be viewed as a correlation function of twist operators on a 
single sheet. 


Tr Paub ~ = (^■( 2 : 1 , Zi)a{z 2 , Z 2 )a{z 3 , Z3)a{z4, 2:4))strip , (2.8) 

where ( 2 ^ 1 ,^!) and ( 2 ^ 2 , ^ 2 ) are the endpoints of region A, and (^ 3 ,^ 3 ) and (^ 4 ,^ 4 ) are 
the endpoints of region B. This correlator is computed on a strip of height (3/2, with 
















boundary condition B at the top and bottom. The twist operators a and d, which 
have opposite orientation, are dehned to reproduce the multisheeted path integral in 
(2.7). This dehnition, together with the conformal transformation properties of the 
path integral, lead to the identihcation of the conformal weights of the twist operators 
(I-o, ro) ^ (hn, /in) as [16,17] 



(2.9) 


The twist operators are not local operators in the original CFT, but they are local 
operators in the orbifold CFT"'/Z„, so we can use standard methods to compute these 
correlators in the orbifold theory. 


For Im Zi = |, with Zi = z*, the correlator (2.8) computes the replica partition 


function at f = 0. To calculate the Renyi entropy as a function of real time, we 
hrst compute the correlation function for Euclidean insertions off the real line with 
arbitrary Zt, then analytically continue to Lorentzian signature by taking z and z to 
be independent complex numbers. For the conhguration (2.1), 


zi = -t + ifi/A, 

Z2 = L - t + 

z^ = D -\- L — f + f/?/4, 

Zi = D + 2L — t + if] / A, 


zi = t — if3/A 
Z2 = L + t — if] / A 
z^ = D L 1 — if]/A 
Zi = D + 2L + t — if] / A 


( 2 . 10 ) 


Mapping to the upper half plane 

The correlation function on a strip (2.8), with boundary condition B at the top and 
bottom, can be conformally mapped to a correlator on the upper half plane (UHP) 
with boundary condition B on the real line, via 


w = e 


.2Trzl^ 


( 2 . 11 ) 


The correlator is then 



( 2 . 12 ) 
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where 


= {a{wi,Wi)a{w2,W2)(T{w:i,Wz)a{wi,Wi))ijBP ■ (2.13) 

In general, is not fixed by conformal invariance and cannot be computed. However 
its properties under conformal transformations can be understood from the method of 
images: The conformal Ward identity for a 4-point function in the UHP is identical 
to the conformal Ward identity for a holomorphic 8-point function in the full plane, 
where for each insertion Wi in the UHP we include an image point 6 , 7,8 = 
in the lower half-plane [18]. That is, '?hi, 2 , 3 , 4 ) behaves under conformal 

transformations exactly like an 8-point function 

(0(Wl)0(M;2)0(M.'3)0(tf.'4)0(w5)0(tf.'6)0(w7)0(tf^8))plane (2.14) 

with the image points 


Ws = ^4 , W(i = W 3 , W7 = W 2 , Ws = Wi. (2.15) 

This is illustrated in figure 2 (for t purely imaginary). We have written 0 rather than a 
in this correlator because the 4-point function on the UHP is not actually equal to an 
8-point function of twist operators on the full plane, it just has the same transformation 
properties. There is not necessarily an actual (local) operator 0 in the theory whose 
8-point function reproduces Juhp. 

2.2 Entanglement entropy in the boundary state model 

The time-dependent entanglement entropy in this model has been discussed many 
times in the literature [6-8,14]. However, one of our goals is to understand exactly 
when those results apply, so we will repeat the calculation in detail. We begin with a 
discussion of the regimes away from the dip, where we agree with the literature, then 
return to the possibility of a dip in the next subsection. 

The reason we can make further progress in the scaling limit t,L,D S> (3 is that, in 
some cases, the correlator in this limit is fixed by the short-distance behavior of the twist 
operators. There are two limits where this short distance behavior is known [8,17,19]. 
First, a twist operator a in the UHP may approach another twist operator a in the 
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Figure 2: Schematic illustration of the method of images. The conformal transfor¬ 
mation properties of the correlator with insertions tci, 2 , 3,4 in the UHP are identical 
to those of the holomorphic 8-point function in the full plane, with image points 
'?^i= 5 , 6 , 7,8 = wg-i. The solid red lines indicate the regions A and B. 


UHP. In this case we can ignore the boundary, and the short-distance behavior is 

a{x,x)a{y,y) {x — {as x ^ y ,x ^ y) . (2.16) 


Second, a twist operator in the upper half plane may approach the boundary. It is 
therefore approaching its image point, and the short distance behavior is 

a{x,x) rn{x — x)~‘^^”' (as x —>■ F) , (2.17) 


where is a constant (since we have already hxed the normalization of the operator 
in (2.16)). 

In Euclidean signature, (2.16) and (2.17) are the only types of short distance limit. 
The singularities in Lorentzian signature are more complicated, since w and w are 
independent; we will return to this below. 

The two Euclidean limits £x the behavior of the Lorentzian correlator in 

certain regimes where all four twist operators are in one of a limits (2.16) or (2.17). To 
see when this occurs, we will look at the cross ratios 


Vij Vij 


{Wj - Wj){Wj - Wj) 
{Wi - Wj){Wj - Wi) 


(2.18) 


in different regimes. Only five of these six real cross ratios are independent, but it is 
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useful to keep track of all of them. 


Early times 

At early times, 


the —>■ 0 limit sends 


L 

0 <t < — , 
2 ’ 

rjij -)■ 1 . 


(2.19) 

( 2 . 20 ) 


This is a conhguration where each twist operator approaches its image point, tCj —)■ tCj. 
This is not immediately obvious, since there are multiple simultaneous limits. For 
example. 


W 2 - Wi 
W 2 - W 2 


—)■ 0 , 


( 2 . 21 ) 


but this does not contradict W 2 —)■ W 2 , as a{wi,wi) goes to the boundary faster than 
to any other operator. That is, wi approaches wi without any other operators nearby 
and this means we can apply the boundary OPE (2.17), a{wi,wi) rn{wi-Wi) 
Having eliminated cr(wi, Wi), we can then repeat the same argument for W 2 —)■ W 2 , then 
for W 3 ^ Ws, and then for w 4 ^. Therefore the leading term in the correlator is 


rUHP 


[(wi - Wi){w2 - W2){W2, - W^){W4 - ^ 4 )] 


— \-\—2hn 


r„exp 


■ ——hn{D + 2L + 2t) 


Plugging into (2.12), 


Jf-P = rtij] exp 


Ibvr 


h„t 


( 2 . 22 ) 


(2.23) 


The resulting entanglement entropy in the scaling limit is 

dvrcf 


c®3-i'h_o Q 

^AUB — ^>^0 + 


3/3 ’ 


where 


2 So = SAusit = 0) = ^ log , 

3 27re 


(2.24) 


(2.25) 


is a constant and we have made explicit the dependence on a UV-cutoff e. 
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Late times 


Now consider 

t>j + L. (2.26) 

In this regime, the /? —?• 0 limit sends 


Vij 0 ) 


(2.27) 


with 


Vl2 ~ h34 < h23 < h24 ~ hl3 < hl4 • 


(2.28) 


In this limit, W 2 —)■ Wi and W 2 —>■ Wi faster than any other limit, so we can apply 
the OPE {w 2 ,W 2 ) —)■ {wi,wi) in the manner of Eq. (2.16). This leaves the operators 
at points 3 and 4, and since 7734 —)■ 0 we can apply the OPE ( 103 , 103 ) —)■ ( 104 ^, 104 ). 
Therefore 


rUHP 


[(Wl - W2){W3 - W4){wi - W2){W3 - ^ 4 )] 

exp -^hn{D + ?,L) 


—2h„ 


(2.29) 


This gives the entanglement entropy 


= 2So + ^. (2.30) 

The second term is the thermal valne, (2L) x Sthermai, where Sthermai is the thermal 
entropy density at inverse temperatnre f3. 

Intermediate times, before the dip 

Now consider 

(2,31) 

The npper limit restricts to the time before the possible dip. In this case, 

rji2 —)• 0, r]i3 —>■ 1, rii4 —1, 7723 —^ 1, V 34 0 . (2.32) 
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Once again we can use the OPEs { 102 , 102 ) —)■ (wi,wi) and then { 104 , 104 } —)■ { 103 , 103 ), 
so the answer is the same as at late times: 


SZb = 25 o + — . ( 2 . 33 ) 

2.3 Light cone singularities and the dip 

Finally we turn to the dip regime, 

f < i < f + i . (2.34) 

In this regime, for t < ^ + f, the limit /? —?• 0 corresponds to 

hi 2 0, r]i3 —>■ 1, r]i4 —)■ 1, rj23 —)■ 0, 7724 —>■ 1, 7734 —)■ 0 , (2.35) 

and for t Y 

hi2 0, 7]i3 —)■ 0, ?7i4 —>■ 1, 7723 0, 7724 —>■ 0, 7734 —)■ 0 . (2.36) 

Now we come to a key point: the conhgurations in eqs. (2.35) or (2.36) do not corre¬ 

spond to any combination of the OPE limits (2.16) and (2.17). Instead, they correspond 
to the limits 


Wi -)■ W2, 
wi -)■ The, 


W3 -)■ Ws, 
W2 -)■ The, 


W4 —)■ Wt, 
W3 -)■ W4, 


W 5 -)■ We 

Wy ^ Ws ■ 


(2.37) 


This is not possible with the OPE, because the Wi are in a different channel than 
the Thj. These limits are therefore intrinsically Lorentzian. They correspond to an 
operator hitting the light cones of two other operators, simultaneously. To illustrate 
this, consider the cross ratio 

(tUi - W2){W5 - We) 

X — - 

{wi - W5){W2 - We) 

^( 2 . 36 ) illustrates why we need to be careful about the order of limits — it is a different limit 
than ( 2 . 27 ), but would seem to be identical if we only looked at the five independent cross ratios 
(7712,7713,7723,7724,7734) and ignored the order of limits. 


(2.38) 
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In the regime (2.34), 

X —0 , X ^ 1 . (2.39) 

If we map to wi = wi = 0,^5 = = 1,wq = = cx), this corresponds to the limit 

where { 102 , 102 ) approaches the tip of the causal diamond bounded by the light cones 
of operators 1 and 5, as in hgure 3. 

In general, the OPE does not £x the behavior of a correlator in this type of limit, 
and in particular it does not imply that there is a singularity [a;(l — On the 

other hand, in rational OFT, these singularities do exist. Much of the rest of the 
paper will be devoted to studying these singularities in detail, and understanding the 
difference in the non-rational case.^ 

To see when we do get a strong dip in the entanglement entropy, suppose that we 
can treat w and w independently, so that we can choose different channels for left and 
right movers. Then the limit (2.37) would produce a term 

~ [{Wi-W2){w3-Ws){w4 -W7 ){w5-Wq) (2.40) 

x(m}i - We){w2 - W5){W3 - W4){W7 - Ws)]"^" • 


This dominates over the other channels in the regime (2.34), and leads to an entangle¬ 
ment entropy 


= 2^0 + g 


D + 2L-2t 
2t-D 


f <^<f+ f 

^+^<t<^+L 


(2.41) 


^ In [6,8] it was assumed that the singularity [a;(l — always exists. In the notation of [8], 

this assumption is implicit in the claim that J^n,N can be set to one in the limit /3 —>■ 0. This is the 
technical origin of the difference between our conclusions and previous analyses. 
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This reproduces the dip of [6-8], plotted in hgure 1. 

To summarize: (i) this dip exists if and only if the 4-point function of twist operators 
in the UHP has the light cone singularity (2.40), and (ii) this singularity is not required 
by the OPE (though we will see that it is required by crossing symmetry in rational 
CFT). 

Note that the additional light cone singularity we are referring to is not simply the 
singularity where an operator in the UHP hits the image point of another operator. 
This type of singularity, which is also not required by OPE, is not sufficient to produce 
the quasiparticle dip. 

3 Light cone singularities in the Renyi entropy 

3.1 Thermal double model 

The method of images determines the conformal transformation properties of a corre¬ 
lator on the UHP, but it cannot be used to actually compute the correlator. In other 
words, the 4-point function of twist operators on the UHP is not fixed by the 8-point 
function with image points of twist operators on the full plane. In particular it depends 
on the specific boundary condition. To avoid this complication, we can use a slightly 
different model for the quench introduced in [5]: simply replace the image points by 
actual twist operators, i.e., compute the 8-point function 

{a{wi,Wi)a{w 2 ,W 2 )cr{w 3 , W3)a{w4, W4)a{w5,W5)a{wG, We)a{w 7 , U)7)d(tC8, u)8))plane ■ (3.1) 

Physically, this calculates the Renyi entropy in a doubled CFT. The doubled CFT 
consists of two decoupled copies of the original theory, CFTi x CFT 2 , in the thermofield 
double entangled state 

|TFD) = . (3.2) 

n 

Each CFT is labeled by its own coordinates, {ti, xi) or (^ 2 , X 2 ). The regions now include 
a piece in each system: 

H = HiUH2, R = RiUR2 (3.3) 

where Hj=i 2 is the interval Xi G [0,T] in system i, and Bi is the interval x* G [D + 
L, D + 2L]. See figure 4a. The CFTs are evolved in time under the total Hamiltonian 
H = H, + H 2 . 
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B 2 



A' 


(a) (b) 

Figure 4: (a) Entanglement regions in the thermal double model, corresponding to 
Saub in the original CFT. (b) Simplihed setup, with an interval in the CFT and an 
offset interval in the thermal double. This setup exhibits the same features, including 
long-term memory effects in the quasiparticle picture that are not necessarily present 
in CFT. 

The Renyi entropies are evaluated at = ^2 = t and have a non-trivial time 
dependence in the |TFD) state. Similarly to what we discussed in sec. 2.1 for the 
boundary state model, they can be obtained from analytic continuation of the twist 
correlator for Euclidean insertions on a thermal cylinder of periodicity j3. The twist 
insertions are at the endpoints of Ai, Bi with a shift 1(3/2 for operators in two different 
copies of the CFT. The 8-point function on the plane (3.1) is related via conformal 
mapping to this 8-point function on the cylinder. 

This thermal double model was designed partly to reproduce the physics of the 
Calabrese-Cardy boundary state model for a single interval. For multiple intervals, it 
is not clear that these models are equivalent. If anything, it appears that (3.1) may 
have more singularities than the corresponding quantity in the boundary state 
model, since we do not know of any argument that the boundary state correlator has 
singularities when insertions in the UHP hit image points of different operators in the 
lower half plane (for example Wi —)■ w^). Therefore we suspect, but have not proved, 
that the Calabrese-Cardy model cannot have any singularities beyond those in the 
thermal double. In this section we will determine what distinguishes the CFTs with 
and without a large dip in the Renyi entropies, in the thermal double model. With 
the (in our view conservative) additional assumption that boundary correlators do not 
have any additional singularities beyond those present in the 8-point function, the same 
conclusions apply to the boundary state model. 

3.2 Offset intervals in the thermal double 

We have argued that entanglement memory comes from the light cone singularity of 
the cross-ratio x introduced in (2.38). To isolate this physics, we will study the simpler 
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correlator 


«’<j)>plane , 


(3.4) 


with the same points as dehned above, 



(3.5) 


If there is a dip in the entanglement compnted from the 8-point fnnction, then it 
mnst also appear in this 4-point fnnction. So, if we can rnle ont qnasiparticle-like 
singnlarities of this 4-point fnnction, then the qnasiparticle pictnre is rnled ont for the 
8-point fnnction as well. 

The correlator (3.4) has a natnral physical interpretation in the thermal donble 
theory: It compntes the Renyi entropy of a region A' consisting of an interval in CFTi, 
and an offset interval in the thermal donble CFT 2 , as in hgnre 46: 


A' : {xi G [0, T]} U {x2 G \D -\- D -\- 2//]} . 


(3.6) 


Namely 



(3.7) 


The qnasiparticle pictnre predicts a memory effect in this sitnation. In the thermal 
donble, the initial state in the qnasiparticle pictnre consists of entangled pairs, where 
one member of each pair is in CFTi and the other is in CFT 2 , with both particles at 
the same x-position. All of these pairs contribnte to the entanglement entropy at t = 0, 
so the initial entanglement entropy is (2L) x Sthermai(/d) (plus the nsnal divergent term). 
Under time evolntion, the particle in CFTi moves one direction, and the particle in 
CFT 2 moves the other direction. In the range 



(3.8) 


we therefore hnd particles in CFTi with xi G [0, T] that are entangled with particles 
in CFT 2 with X2 G [D + L,D + 2V\. This is illnstrated in hgnre 5a. Both of these 
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(a) (b) 

Figure 5: Quasiparticle picture for the offset-intervals entanglement in the thermal 
double model, (a) Entangled quasiparticle pairs (blue dots) in systems 1 and 2 are 
separated by Ax ~ /3 at t = 0 (dashed horizontal line). Under time evolution, the 
particles move in opposite directions. Eventually both particles enter region A', (b) 
Quasiparticle prediction for the entanglement entropy Sa' — 2S'o of region A'. When 
the entangled pair enters region W, the entanglement entropy decreases. 


particles are in region W, so the entanglement entropy decreases. Thus the quasiparticle 
prediction for the entanglement entropy, corresponding to the u —>■ 1 limit of (3.4), is 
as shown in hgure 56. The early-time prediction is exactly computed for 

the boundary state in (2.30), and the dip prediction is the same as in (2.41). Our 

goal in this section is to understand when this quasiparticle prediction is correct and 

when it is incorrect. 

The correlator (3.4) is (for f, L, H S> /3) 

{xxY’^"Gn{x, x) , (3.9) 

where 

Gn{x,x) = {a{0)a{x, x)a{l)a{oo)) . (3.10) 

The cross-ratio x, dehned in (2.38), is 


^plane _ 


exp 


■jK{D + 3L) 


X 

X 


exp 

exp 




2ti 


— {D + 2L + t — max(Zl -|- 2L 


t,t) 


max(Zl, 2t)) 


(3.11) 
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The OPE implies that this correlator has the singularities 


Gn ~ (xx) (as x, a; —)■ 0) 


(3.12) 


and 

~ [(1 — x)(l — (asx,x—)-l). (3.13) 

Away from the dip, i.e., outside the range (3.8), x ~ 0 so the correlator is given by 
(3.12). This leads to the entropy 

yearly, late ^ ^ ^ (3.14) 

as expected from the quasiparticle picture. 

In the range (3.8), x —0 and x —)■ 1, with 


1 — X ~ exp 


2vr 

■— min(Zl + 2L 

P 


2t, 2t - D) 


(3.15) 


so that 

0 . (3.16) 

1 — X 

The OPE does not require any particular behavior in this limit. But if we assume that 
Gn{x, x) ~ [x(l —x)]“^^" in this limit, then the resulting entanglement entropy is equal 
to (2.41). This would agree with the quasiparticle picture. Therefore in this simple 
setup, the validity of the quasiparticle model is precisely the question of whether the 
correlator (3.10) has a singularity 


G„(x, x) ~ [x(l — x)] (as X —)■ 0, X —>■ 1). (3.17) 


3.3 Light cone features in the second Renyi 

Our conclusion will be that the singularity (3.17) exists if and only if the theory is 
current-dominated, as explained in the introduction. It is simplest to make this argu¬ 
ment for the second Renyi entropy, 

G2 (x,x) = (cr2(0)d2(x,x)cr2(l)d2(oo)) , (3.18) 
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where the subscript indicates replica number n = 2. The two-sheeted manifold that 
dehnes these double-twist operators is in fact a torus. It follows that this correlation 
function is hxed by the torus partition function, and the mapping derived in [19,20] is 

G 2 {x,x) = {2}^xx{l — x){l — x)) (3.19) 


The partition function. 


Z(r,f) = Tr exp 


2'KiT ( ^0 - - 27rir [Lq - ^ 


(3.20) 


is evaluated with the torus modulus related to the twist operator insertions by 


X = 


Mil! 


T = l- 


K{l-x) 


K{x) 


(3.21) 


where K{x) = | l,a;^), and similarly for t{x). For small x, the relation is 


X = I6y/q + 0{q), q = e 


_ Jl'Kir 


(3.22) 


and for x near 1, 

1-a; = 16^ + 0(g'), g'= . (3.23) 

The modular transformation of the torus r —)■ — 1/r is related to crossing symmetry of 
the 4-point function, x ^ 1 — x. 

This maps the problem of memory effects in the 2nd Renyi to a question about the 
spectrum. Including the prefactors from (3.19) and (3.9), and the conformal factor in 
(3.7) from mapping the cylinder to the plane, the 2nd Renyi is 

Sj = ^ - log [2-^l\xxrl^\(l-x)[l-x)]-^I^^Z(T(x)Mx))] + . (3.24) 

where 

is a constant independent of t^L,D. 

First, let’s reproduce the OPE singularities (3.12) and (3.13). The limit x,x ^ Q 
is the zero temperature limit of the partition function r —>■ ioo,f —)■ —ioo, so only the 
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vacuum contributes to the trace: 


Z{t, f) ~ _ 2^cIZ^-cI12--cI12 _ 

Including the prefactor in (3.19), 

G 2 ~ . (3.27) 

The other OPE limit, a:, a: —?• 1, is a high temperature limit. In this limit, we use 
modular invariance of the torus partition function, projecting the trace onto the vacuum 
in the dual channel 

Z{T,f) ^ (3.28) 

which leads to 

G2~[(1-x)(1-7)]-2'^^ , (3.29) 

as expected. 

Now consider the mixed limit a: 1 — a: 1. This is a mixed-temperature limit, 

where we send 

r —)■ ioo, f —>■ i0~ . (3.30) 

This is low temperature for left movers, and high temperature for right movers. We 
first take the limit a: —?• 0 , so that the partition function gets contributions only from 
left-moving ground states, i.e., conserved currents, 

Z(r,f) ~ . (3.31) 

currents 

The sum is over states with conformal weights (0,Lo). The limit that is actually 
relevant to the 2nd Renyi simultaneously takes x —?• 0 and x —)■ 1 in a particular ratio, 
so (3.31) might not contain the dominant terms in that limit. However our goal is only 
to check whether the quasiparticle singularity is present, and (3.31) is sufficient for this 
purpose. 

Now in the limit x —>■ 1, the sum is dominated by the states with very large Lq. The 
density of states in CFT, Lq), is hxed asymptotically by the Cardy formula [21], 

H(Lo,To) ?»exp f27rW^Lo + 27rW^Loj (as Lq, To 00 ) , (3.32) 
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where c is the central charge. This formula applies only when both Lq and Lq are large, 
so it does not £x the behavior of (3.31). Instead, define the density of currents 

^currents (-^o) = ^(-^0 = 0)-^o) • (3.33) 


The singularity in (3.31) will be fixed by the asymptotic growth of flcurrents- Let us 
parameterize this growth by a number Ccurrents, defined such that 

f^currents (Lq) ~ CXp 

Clearly c > Ccurrents- To see that we can always parameterize the asymptotic growth of 
f^currents in this Way, hrst note that the growth must be at least this fast, since the stress 
tensor guarantees Ccurrents ^ Also, the exponent cannot have a larger power of Lq, 
since a singularity stronger than the quasiparticle singularity would not be consistent 
with modular invariance. 

Ccurrents cau be interpreted as an effective central charge for the conserved right- 
moving current sector of the theory. As g —)■ 1, the sum (3.31) can be evaluated by a 
saddlepoint approximation, with the usual Cardy result 


n / Ccurrents"^ 

^ / ---To 


as 


Lq —)■ oo 


(3.34) 


^ 7") Q ^^^^(^0 2(^^^f^rirrents)/3^ Ccurrents/12 


(3.35) 


This formula is one of our main results: it implies that the quasiparticle picture for 
the second Renyi applies if and only if Ccurrents = c. That is, the quasiparticle picture 
holds only in theories which are current dominated, meaning the asymptotic number 
of conserved currents is approximately equal to the total number of states. This is true 
in rational CFT (see below), but not true in general. In particular, a theory with no 
conserved currents besides the stress tensor has 


Ccurrents A 1 


(3.36) 


Thus, in a theory with c > 1 and no extended chiral symmetry, the quasiparticle picture 

^Without null states, the stress tensor descendants of the vacuum grow asymptotically as ^ 
exp(7r-\/2Lo/3), corresponding to Ccurrents = 1- With null states this number can decrease, with the 
lower bound Ccurrents = 5 Set by the Ising model since this is the unitary minimal model of lowest 
central charge. 
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does not produce the correct 2nd Renyi entropy."^ 

In a non-rational theory with Ccurrents < c, the 2nd Renyi entropy is not universal, 
since representations other than the vacuum may dominate the partition function. The 
contribution from the vacuum representation, obtained by inserting (3.35) into (3.24), 
sets an upper bound during the dip regime, and the quasiparticle result sets a lower 
bound: 

^ ^ min(T) + 2L - 2t, 2f - D) + (3.37) 

< ^ ^ + mm{D + 2L- 2t, 2t - D) + . 

For Ccurrents = c, the dip in S 2 is as large as possible; for Ccurrents -C c there is still a dip 
but its magnitude is smaller by a factor of at most 3. 

More comments on rational vs general CFTs 

In a rational CFT such as the c < 1 minimal models, the partition function is a hnite 
sum of characters:^ 

^(Ri") = '^Xr{(1)Xr{q) ■ (3.38) 

R 

These characters transform under S':r—)■—1/rby the action of the modular R-matrix, 

Xr{q) = SRPXpjq) . (3.39) 

p 

This is again a hnite sum. It follows that the mixed-temperature limit of the partition 
function relevant to the Renyi with a:-Cl — x<^lis 

Z(r,f)~5oog-'^/''(r)"^/'". (3.40) 

Comparing to (3.35), we see c = Ccurrents in rational CFT. 

In a c > 1 CFT where the vacuum representation has only Virasoro descen¬ 
dants, the manner in which modular invariance is enforced is quite different. The 

^In the Wn minimal models, which generalize the Virasoro minimal models to a Wn chiral algebra, 
Ccurrents < IV — 1. In a theory with chiral algebra Wn, but c > iV — 1, there are no null states in the 
vacuum representation and this leads to Ccurrents = iV — 1. 

®We assume the diagonal modular invariant but the conclusions are similar for other invariants. 
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S'-transformation of the vacuum character is an integral [22]® 


f 

xTm = / 

J (c 


'OO 


dhS(h,0)xr(q') . 


( 3 . 41 ) 


(c-l)/24 


This is clearly not a sum over characters in the dual channel, and does not include a 
contribution from the vacuum. As g —)■ 1, 



(3.42) 


so Ccurrents = 1- Modular invariaucc still holds, but in a non-rational CFT, the vacuum 
singularity as r —>■ z0“ does not come from any individual character in the original 
trace. Instead it comes from the asymptotics of the inhnite sum. 

3.4 Conformal block expansion 

We now return to the general Renyi index n. In this case we can make a similar 
argument using conformal blocks instead of characters, but the singularity in the dual 
channel is poorly understood. 

The twist correlator may be expanded in conformal blocks in any channel. In the 
s-channel a: —?• 0, 



(3.43) 


p 


where Op is a constant related to the OPE coefficients, and Fp{d, h] x) is the conformal 
block with internal weight hp, external weights h, and central charge c! = cn (since this 
is a correlator in the replica theory). The conformal blocks are, by dehnition, hxed 
entirely by the chiral algebra. In the limit a; 1 — x -C 1, the dominant contribution 
is not necessarily the vacuum. The contribution from the vacuum sets a lower bound. 



This contribution is universal, in the sense that it does not depend on the full details 
of the theory: it depends only on the conformal block Fo{cn, hn] x) in the limit x —>■ 1. 


6 


S{h,0) 



sinh(27r6 ^ph) where h = p\ + 


c-l 

24 


and c=l+ 6(6 + 5 ^)^. 
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It is therefore fixed by the right-moving chiral algebra. Of course, if the two intervals 
are offset in the opposite direction, then left-moving currents contribute instead. In 
the original quench model, with two intervals in a pure state of a single CFT, both 
types of currents contribute. 

The conformal block expansion (3.43) is in the orbifold, CFT”/Z„, so Fq in (3.44) 
is the vacuum block in the orbifold. The chiral algebra of the orbifold is hxed by 
the chiral algebra of original CFT, but is not identical (see for example [19]). This 
distinction is discussed further below. 

Rational CFT has the quasiparticle dip 

In a rational CFT, we can expand any individual conformal block as a hnite sum over 
primaries in the dual channel; 


~ ’ ( 3 . 45 ) 

g 

(with some of the arguments of the conformal block suppressed). Every term con¬ 
tributes on the right hand side. Thus 

Fq{x) ~ 6oo(l “ 5;)“^^" as h —>■ 1 . (3.46) 

This is the quasiparticle singularity (3.17).^ Thus in the scaling limit f, L, H the 
quasiparticle picture does apply to rational CFTs. 

One way to understand the origin of the singularity in Fq[x) as x —>■ 1 in rational 
CFT is from crossing symmetry, 

Gn{x,x) = Gn{l - X,l - x) . (3.47) 

In a rational CFT, each side of this equation is a hnite sum. Therefore, in order to 
reproduce the identity singularity on the rhs, individual terms on the Ihs must go as 
(1 — 

^The coefficient &oo in (3.46) is related to the quantum dimension of the twist operator in the 
orbifold theory. It does not affect the answer in the scaling limit t,L,D ^ f3, but may have inter¬ 
esting implications at intermediate times; quantum dimensions have appeared several times before in 
entanglement calculations [23-25]. 
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General CFT 


In a general CFT, the conformal block expansion is infinite. In this case the crossing 
eqnation (3.47) does not imply anything abont the singularities of individual terms 
in the s-channel as x —)■ 1, since these singularities can come from the inhnite sum 
rather than from any particular term. This is exactly what occurs in higher spacetime 
dimensions, where there is no such thing as rational CFT — individual conformal 
blocks cannot reproduce the identity singularity in the other channel, so it is instead 
reproduced by the asymptotics of the inhnite sum (see for example [26]). 

Suppose that the chiral algebra of the CFT is just the Virasoro algebra. That is, 
the only states in the CFT of dimension (0, Lq) are the Virasoro descendants of the 
vacuum. The vacuum contribution to the twist correlator for x-Cl — x-Clis then 

Gn{x,x) ~ ^^"(cn,/i„;x) , (3.48) 

where is the vacuum conformal block associated to the chiral algebra N'vF j^Ln- 

This conformal block is not just the Virasoro block, though it is hxed entirely by the 
Virasoro algebra. Denoting the Virasoro modes in copy k of the CFT hj \ the block 
^Vir /z„ dehned to include the contributions from all Z„-symmetric combinations of 
the states 

( 3 . 49 ) 

To proceed, we would need the behavior of this vacuum block in the limit x —)■ 1, but 
this is not known. Indeed, even the behavior of in the limit x —)■ 1 appears to be 
unknown. However, there is generically no reason for it to reproduce the quasiparticle 
singularity (3.17). We therefore expect in this limit 

Fo^""/^"(cn, hn,x) < (1 - x)-"'*" . (3.50) 

Our results for the second Renyi entropy in section 3.3 conhrmed (3.50) for n = 2. We 
do not have a conclusive argument argument for n > 2, but since the singularity is 
not required it is highly implausible that it would appear (for generic central charge 
— of course it does appear in the minimal models). Crossing symmetry does not allow 
a singularity stronger than (1 — x)“^^", so this leaves (3.50). Other than the second 
Renyi, there is one other case we are aware of with c > 1 where the Virasoro blocks can 
be computed exactly (c = 25, h = 15/16 [27]), and the inequality analogous to (3.50) 
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holds in that case as well.® 


Since we do not know the singularity of the orbifold conformal block as a function 
of n, we cannot determine the entanglement entropy, beyond the statement that it 
disagrees with the quasiparticle picture. It would be very interesting to compute the 
universal vacuum contribution. 


3.5 c=l 

We have discussed the minimal models with c < 1, and “generic” CFTs with c > 1 
but no conserved currents other than the stress tensor. Now we turn to theories with 
c = 1 and show that the Renyi entropies of these theories agree with the quasiparticle 
picture. 

The starting point for almost all known c = 1 CFTs is the compact free boson 
theory,® dehned by the Lagrangian 

L = ^dXBX , (3.53) 


with the compactihcation condition X = X + 27rR, i.e., the target space is a circle with 
radius R. This theory has c = 1 and infinitely many (Virasoro) conformal primaries 
with a discrete spectrum of conformal weights for hnite R [29,30]. 

We can analyze the light cone singularity of the second Reyni entropy of this theory 
explicitly, via the torus partition function [30] 


Z{R) 


1 

ri{r)fj{f) 





(3.54) 


^Explicitly: 


Fo ( c = 25,/i= ) = 


lQq{x{l - 


1 -1 


q = exp 


_K{l-x) 

K{x) 


(3.51) 


where 63 and K are the standard elliptic functions. In the limit x —)■ 1, we have q ^ e i°g(i6/(i x)) 
and 03 (g) ~ V ignoring constants and log corrections, 


Fn ( c = 25, /i = ^ 


; x^ ^ (1 — x) (as X —>■ 1) . 


(3.52) 


This is much weaker than the vacuum singularity (1 — x) = (1 — x) 
®The theory in [28] is an exception. 









where the conformal weights are given by 


1/e mR\^ 
^ 2\R^ ^) 

1/e mR\^ 

^ 2[r~ ^) 


(3.55) 

(3.56) 


and r] denotes the Dedekind eta function. When R^ is rational, there are inhnitely 
conformal primaries with a given he^m or he,m- For example, = 0 whenever 2e/m = 
—R^. Moreover, in this case the Virasoro algebra is extended to a larger chiral algebra, 
with respect to which there are hnitely many primary operators, making the theory 
rational [30,31]. But when R^ is irrational, there is no such degeneracy in the weights 
and the theory has no extended symmetry, so we investigate this case separately (these 
CFTs are examples of what [32] calls “quasi-rational”). 

In the light cone limit r —)■ too, f i0~, q —>■ 0+, g —)■ 1“, the partition function in 

Eq. (3.54) behaves as 

Z{R) = , H- ] , (3.57) 

v{T)v{r) V ) 


where hmm and are the minimal non-zero holomorphic and anti-holomorphic 
conformal weights, respectively. The assumption that R^ is irrational ensures that 
hmin, hjnin > 0, 1 . 6 ., there are no additional conserved currents. Therefore, all the terms 
besides 1 approach 0 in this limit, and the series is absolutely convergent, so the be¬ 
havior of the singularity is controlled by the rj functions. To analyze them we use the 
modular transformation property: 


r]{-l/T) = , 



(3.58) 


where r' = —l/r. In this limit r' —)■ —ioo and so we have t]{t) —and rfir') 
and thus 

Z{R) ~ g-i/24^/-i/24 ^ 

where the last expression gives the asymptotic behavior because the q and q divergence 
is exponential and so dominates over the vanishing of f. We thus expect the dip in the 
second Reyni entropy after a quench, because we have the same asymptotic behavior 
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as the rational CFT case, shown in eq. (3.40). 

The Reyni entropies of two intervals in the compact boson theory were computed 
in [33,34], and generalized to arbitrary complex cross ratio in [35]. We will consider 
only the decompactified limit, i? —>■ oo. In this case the correlator of four n-twist 
operators is 


{a{wi,wi)a{w2, ta2)cr(tC3, w^)a{wi, w^)) 


(3.60) 



with X = and Wij = Wi—Wj. The explicit formula for J^n(,x, x) is given in equation 

(142) of [35], from which it is easy to check that the only power-law singularity in the 
light cone limit tC 2 —t tci, tC 2 —t Ws, which corresponds to x —>■ 0, x —>■ 1, comes from 
the prefactor in (3.60). This is the quasiparticle singularity. 

The general picture, then, is that for a given chiral algebra, there is a critical central 
charge Ccurrents- In many examples, including Virasoro symmetry, Wn symmetry, and 
Kac-Moody symmetry (and perhaps in general), it is the same as the threshold for the 
theory to be rational. If c = Ccurrents, including irrational but ‘nearly rational’ theo¬ 
ries like the irrational boson where the central charge sits at the threshold, then the 
spectrum is current dominated and the entanglement entropy has quasiparticle behav¬ 
ior; and if c > Ccurrents then the entanglement entropy disagrees with the quasiparticle 
picture. 

3.6 Entanglement is larger than quasiparticle at the dip mid¬ 


point 


We have argued that in theories with c > Ccurrents, the entanglement entropy of two 
intervals does not agree with the quasiparticle prediction. There is one hnal step before 
we can conclude that memory effects are reduced: we need to argue that the dip is 
smaller in these theories, not larger. This seems obvious physically, but we will only 


give a proof limited to the second half of the dip regime: < t < ^ + L. This is 

enough to demonstrate the result at the midpoint of the quasiparticle dip: 




(3.61) 
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To demonstrate this, we apply strong subadditivity to the neighboring regions 


A:[0,L], C :[L + 5,L + D-5], B : [L + D,2L + D] , (3.62) 

separated by a small distance, S /3 t, L, D. The entanglement entropy of each 
individual region is given by the Calabrese-Cardy formula for a single interval of size i 
after a global quench: 

TTC 

Se{t) = 5*0 + — min(2t,£) . (3.63) 

We can also apply their multiple-interval result to neighboring regions (since for separa¬ 
tions of order 6 contributions from possible light cone singularities would be suppressed 
as 5/P ): 

TTC 

5'auc = ScuB = 2*^0 + min(2t, L + D) , (3.64) 

and 

TTC 

Saubuc ^ ^ min(2f, 2L + D) . (3.65) 

Now strong subadditivity, in the form 

Saub > ‘S'^ucuB + Sb — ScuB (3.66) 


implies 


TTC 


S'aub > 2S'o-I—- [min(2f, 2L-|-D)-|-min(2f, L) — min(2f, L-|-D)] . (3.67) 

ojD 


In the range < t < ^ + L, this becomes 


Saub > 2S„ + ^(2« - D) = 53“”“““ ■ (3-68) 

Thus the quasiparticle prediction is a theoretical lower bound for the entanglement 
entropy during the second half of the dip (the rising edge of the dip in hgure 1), in 
particular at the midpoint (3.61). 

It is interesting to note that a similar bound for the hrst half of the dip regime 
f <t<^ could be obtained from the monogamy inequality 


Saub > Saucub + Sa + Sb + Sc — 5'cuu4 — Scub 


(3.69) 
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which is believed to hold in holographic theories [36,37], but is not satished in general 
quantum systems. 

4 The large-c limit 

We have argued that entanglement scrambles for c above a critical value, but the actual 
value of the entanglement entropy in theories that scramble is not universal during the 
dip regime. In this section, we will show that in the holographic c —>■ oo limit, all 
theories scramble maximally. We will derive this universal limit using semiclassical 
gravity, and conhrm it in conformal held theory (subject to some assumptions about 
the growth of OPE coefficients in these theories). 

Holographic CFTs have few chiral states, so Ccurrents ~ 1, but many total states, so 
c 3> 1. Therefore some of the GET calculations in sections 2 and 3 must be revisited, 
since the large-c limit may compete with the scaling limit t,L,D ^ j3. In section 4.4, 
we address the interplay between these two limits (see also [38]) using the example of 
the 2nd Renyi in the D1-D5 GET. 

4.1 Holographic entanglement entropy 

The holographic dual of the thermoheld double setup described in section 3 has been 
studied in [5]. It consists of the eternal BTZ black string with time taken to run 
forwards on both sides of the Penrose diagram, thus producing a time dependent excited 
state. For convenience, in this section we will set the inverse temperature {3 of the black 
string to 27r, and only reinstate (3 in the hnal expression of the holographic entanglement 
entropy. We also set the AdS radius to ^Ads = 1- 

According to the proposal of [39,40], the holographic entanglement entropy of dis¬ 
joint intervals in 2d GFTs is computed by the minimal length collection of geodesics 
that extend in the bulk and join the intervals endpoints on the boundary [19]. To com¬ 
pute the geodesics associated to the offset intervals, we follow [5] and exploit the local 
equivalence between the BTZ black hole and empty AdSs. This allows to compute the 
length of geodesics in the black string geometry directly in terms of the length of AdSs 
geodesics anchored on the conformal boundary, which are just semicircles. 

As reviewed in [5], the different regions of the BTZ black string can be mapped to 
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the Poincare patch of AdSa: 


ds^ = \{-dyl + dyl + du^). (4.1) 

In particular, the exterior metric of the black string 

= — sinh^ p dtl + cosh^ p dx\ + dp^ , (4.2) 

with conformal boundary at p —)■ oo, maps to a portion of the Poincare patch (we refer 
to the Appendix B in [5] for the full explicit coordinates transformation). Near the 
boundary, the two sets of coordinates are related by 

Pi ± Po ~ ~ . (4.3) 

u 2 

These coordinates cover one of the exterior regions of the eternal BTZ string; the 
second exterior region with spatial coordinate X 2 is reached by continuing f —)■ f + ztt. 
The boundaries of the two exterior regions thus cover the portion Pq — y\ < ^ oi the 
AdS boundary. In particular, a point P on the right (left) boundary of the black sting 
in Poincare coordinates reads 


p = {yo,yi) = (sinh± cosh fi^a). 


(4.4) 


Now, recall that the length of a spacelike geodesic in AdSs joining two points Pi, 
P 2 at radial coordinates up^ near the AdS boundary m = 0 is given by 


Ci 2 = log 


- Apg + Ayl 

Up^^Up^ 


( 4 . 5 ) 


where Aj/q = ( 2 / 0 ,Pi ~ 2/o,P 2) similarly for Ayi. Using (4.3), it is then straightforward 
to obtain the length of spacelike geodesics anchored on the boundaries of the black 
string geometry. 

The offset interval setup we are interested in corresponds to one interval of length 
L on each boundary, shifted by a distance D, at some hxed time t. In Poincare 
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Figure 6 : Illustration of the two competing configurations of geodesics for the offset 
interval in the {yi-,u) plane. The disconnected conhguration in solid lines is the one of 
minimal length for all times. 

coordinates, this is the set of pointsd® 


Pi = (sinh f, cosh t) , P 2 = e^(sinht, coshf) , 

P 5 = — coshf), Pq = — coshf) . 


(4.6) 


In the computation of the holographic entanglement entropy, there are two competing 
sets of bulk geodesics ending on the endpoints of the offset interval. One consists of 
disconnected geodesics joining Pi,P 2 and P^^Pq] while the second connected conhgu¬ 
ration of geodesics extends throughout the interior region of the BTZ string joining 
Fi, Fe and F 2 , F 5 (see hg. 6 ). 

To regularize the length of these geodesics we introduce a radial cutoff Pcutos in 
the coordinates (4.2) and dehne e = In Poincare coordinates, the cutoff 

is Mcutofri ,2 = and depends on the spatial coordinate xi ^2 of the point on the 

boundary. From (4.5) and (4.6), we then immediately obtain the entanglement entropy 
associated to the disconnected conhguration 




C 


12 


disconnected 


F 56 2 c 


2 sinh(L/ 2 ) 


(4.7) 


where we have used c = ?)/{2Gn)- Reinserting the correct factors of f3 and taking the 
We use the same notation for the points as in sec. 3.2. 
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scaling limit L ^ /3, we have 




disconnected 


— 250 + 


27rcL 
3/3 ’ 


(4.8) 


where, as in (2.25), we have isolated the divergent contribution 


2So=?^log^, 


(4,9) 


For the second, connected, conhguration we have instead 


q _ _ C , 

*^connected ~A~py o 

4(jr7V O 


cosh(2t) + cosh(Zi) + L) 


(4.10) 


Reinserting the correct factors of (3 and taking t^D^L 3> /3, this reduces to 


5, 


connected 


25o + 


27rc 

3^ 


D + L for P<t<^ 

2 t for t > 


(4.11) 


Minimising over the different conhgurations then shows that, as expected, the discon¬ 
nected conhguration is always dominant, so the entanglement entropy is constant: 


5 


5, 


disconnected 


— 25o -l- 


2ttcL 

3/3 


(4.12) 


for all times. This corresponds to the ‘maximal scrambling’ (dashed) line in hgure 1, 
and to an identically vanishing mutual information. 

In particular, this result cannot be corrected by complex geodesics connecting real 
boundary endpoints, as in pure AdS these have the same length (4.5) of real ones. 


4.2 Second Renyi entropy from semiclassical gravity 

In section 3.3 we studied the second Renyi entropy of the offset interval through the 
analysis of the torus partition function Z{t, f), and argued that generically the quasi¬ 
particle argument does not produce the correct 2nd Renyi. Here we compute the torus 
partition function from semiclassical gravity and explicitly verify that in holographic 
CFTs the quasiparticle dip in the second Renyi of the offset interval is suppressed by 
a factor 3. 
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On the gravity side, the torus partition function is computed by a path integral 
over Euclidean three-geometries that are conformal at inhnity to a torus of modular 
parameter t = ^. The inverse temperature and angular momentum are respectively 
the real and imaginary parts of (3. 

In the semiclassical limit, all known saddlepoints of pure 3d gravity that contribute 
to the partition function have classical action of the form [41,42] 

ZTT 

S{T,f) = (4.13) 

oCjttv 

^eSL{2,Z). (4.14) 

Note that the same S'L(2,Z) transformation acts on both r and f. The simplest 
saddlepoint is thermal AdSs, with action 

ZTT 

^thermal (r, f) = ^^^(r - f) , (4.15) 

which can be viewed as a solid torus with a contractible spatial cycle. Its S'-transformation 
r —- is the Euclidean BTZ black hole 

r 

The more general manifolds that can be constructed by a transformation 7 as in 
eq. (4.13)-(4.14) form an S'L(2,Z) family of black holes. 

The leading semiclassical approximation to the partition function is given by the 
solution of least action: 

\ogZ{T,f) ^ -Smm{T,f) = —^max.^ (-*7 -T + ij-f). (4.17) 

o(jr TV 

This leads to a rich phase diagram in the upper-half r plane. Thermal AdS space is 
the dominant classical solution in the fundamental domain |r| > 1, |Rer| < 1/2, and 
in any of its translates by r —)■ r -|- n for integer n, as the free energy is invariant under 
the T-transformation r —)■ r -|- 1. Similarly, the Euclidean BTZ black hole dominates 
whenever there exists an integer n such that —^ lies in the fundamental domain. 


with 


7.7- = 


ar 


CT + df 


7 = 
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The phase diagram resulting from (4.17) is a subtesselation of the usual tesselation of 
the upper half r-plane by fundamental domains of S'L(2,Z), with an inhnite number 
of phases [42], For instance, for purely real (3 one recovers the familiar Hawking- 
Page transition between thermal AdS at low temperatures and Euclidean BTZ at high 
temperatures. 

The modulus of the torus r is related to the cross ratio x in the four-point function 
that computes the second Renyi entropy of the offset interval through eq. (3.21) 


t{x) 


.K{l-x) 
* K{x) 


(4.18) 


The image of the x plane under the map x —)■ t{x) is shown in hgure 7a. Each region 
on the r plane is labeled by the element of SL{2, Z) that takes this region to the 
fundamental domain. This element also provides the value of the action in that region 
— for example, in the region labeled STS, the action is proportional to the imaginary 
part of STS ■ r = Of all the regions shown, there are only four distinct values 
of the action, so these correspond to four different phases, shown in different colors. 
The Euclidean action relevant to the twist correlator is the value of the action on the 
dominant saddle: 


\ogZ[x, x) 


vr 


80 


N 


I I 


max < —IT S IT, - 3 , 


i i 

+ 


IT IT 

+ 


r r’ 1 — r 1 — r’ 1 + t l-|-rj’ 

(4.19) 

where r, r are functions of x, x. The phase regions on the x plane are shown in hgure 
7b. 

In Lorentzian signature the cross ratio varies inside the interval [0,1]. This comes, 
upon analytic continuation, from the Euclidean diamond 


Rex > |Imx| , 1 — Rex > |Imx 


(4.20) 


in the x—plane, shown as the dashed diamond in £g. 7b. As this Euclidean diamond 
intersects only the 1 and S phases, there is only one possible phase transition in 
Lorentzian signature: 


TIC 


log Z (x, x) ~ — max —ix -1- ix , 



( 4 . 21 ) 
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(a) (b) 

Figure 7: (a) Phases of the partition function on the r plane. Regions are labeled by 
the element of SL{2,'L) which takes points in the region to the fundamental domain, 
marked ‘1’. Distinct phases are the four colored regions, separated by bold lines. (6) 
Phases of Z(r(a:), f (a;)), on the x-plane. Only the distinct phases 1, S, ST, STS are 
labeled. The dashed diamond is responsible for the two phases in Lorentzian signature. 




At early and late times {t < D/2,t > D/2 + L): x,x ^ 0, the 1 phase dominates and 

log Z{x, x) ^ ^4 22) 

in complete agreement with (3.26). Substituting in the expression for the second Renyi 
(3.24), we obtain 

yearly,late ^ ^ ^ gO _ ^^ ^3) 

Ip 

At intermediate times {D/2 < t < D/2 + L), we are in the mixed limit where the 1 
phase dominates again asx-Cl — x<^l. This gives 

\ogZ{x,x) ~ (4.24) 

Ga ~ (4.25) 

and 

^ {D + 2L- 2t, 2t-D) + 5° (4.26) 

in agreement with (3.37) for Ccurrents "C c. (We have dropped a term ~ clog2 which is 
subleading in the scaling limit /3 —)■ 0.) 

Observe in particular that a light cone singularity that leads to a second Renyi 
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in full agreement with the quasiparticle picture could only come from saddlepoints 
with action proportional to r — 4. These are the kind of holomorphically factorizing 
saddles that were sought, but not found, in [42], It was pointed out in [42] that such 
saddlepoints would never dominate the classical action in Euclidean signature, but they 
can dominate in Lorentzian signature. 

4.3 Large-c CFT 

We will now discuss how these results are reproduced in large-c CFT. 

The universality class of CFTs with a semiclassical holographic dual is not charac¬ 
terized precisely, but it must satisfy at least two simple criteria (discussed for example 
in [43,44]). First, these theories have c ^ 1. The central charge is related to gravita¬ 
tional parameters by c = 3£Ads/(2GAr), where ^Ads is the anti-de Sitter radius and 
is the three-dimensional Newton’s constant, so the semiclassical limit on the gravity 
side is the c —>■ cx) limit in CFT. Second, the spectrum of low-dimension operators must 
be sparse — each low-dimension operator corresponds to a light field on the gravity 
side, and the number of such fields should be finite as c —)■ oo. Exactly how to define 
the sparseness condition is not known (except for the partition function [44]), but we 
will assume that the spectrum is sufficiently sparse to suppress certain contributions 
to the conformal block expansion. 

In section 4.1 we computed the holographic entanglement entropy of the offset 
interval in the thermal double model. As discussed in 3.2, the entanglement and Renyi 
entropies are obtained in terms of the four-point twist correlator in the plane Gn{x,x) 
in (3.10) via 

^n = (4.27) 

1 — n 

with 

( R 'TT T \ 

^sinhyj \xx\^^-Gn{x,x). (4.28) 

The general semiclassical conformal block analysis to evaluate Gn(x,x) for CFTs with 
an holographic dual has been performed in [45]. There it was argued that, at leading 
order in c, the four-twist correlator is given by the dominant contribution between the 
s- and t-channel semiclassical identity contributions. That is, in these theories, the 
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Renyi entropies are universal, with the twist correlator given by 


Gn{x,x) = max {Fo^cn, hn] x)Fo{cn, hn] x), Fo(cn, h„;l - x)Fo{cn,hn;l - x)) , 

(4.29) 

with Fq the vacuum block for the algebra ViG/'Ln- For n = 2, a complete derivation 
of this formula with a suitable sparseness assumption is given in [44], and the results 
match exactly with the gravity analysis of the 2nd Renyi in section 4.2. For general 
n, additional caveats apply to (4.29), as discussed in [45]. Assuming that it is correct, 
the right-hand side can be evaluated in the limit n —>■ 1 following [45], with the result 
that only the identity operator contributes: 

Gn{x,x) ^ [minima;, (1 — a;)(l — :r)}]“6*'"'“F _ (4.30) 


There is no quasiparticle singularity. 

Applying this result to the computation of the entanglement entropy of the offset 
intervals, one immediately recovers the analysis of the previous section and the holo¬ 
graphic result (4.12). For all times, in the scaling limit L 3> /3, Gn{x, x) is given by the 
s-channel semiclassical identity block, resulting in the constant entanglement entropy 


S 


2So + 


2ticL 

313 


( 4 . 31 ) 


A similar analysis applies to the 8 -point function in the thermal double model for 
the quench. Once again the n —)• 1 limit of the twist correlator in a large-c CFT is 
given by maximizing the identity contribution over OPE channels. Each OPE channel 
corresponds to a choice of Ryu-Takayanagi geodesics connecting the 8 twist insertions, 
so this CFT result is identical to the gravity prediction. 


4.4 The D1-D5 CFT 

A specihc realization of AdS 3 /CFT 2 duality, and hence an in-principle non-perturbative 
formulation of quantum gravity, is afforded by the D1-D5 CFT. We will use this 
example to illustrate how memory effects in the entanglement depend on N, f3, and 
the CFT coupling constant. 

This system is realized in type IIB string theory compactihed on a circle times 
with either a four-torus or a K3 surface, with Ni Dl-branes wrapped around 
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the circle and D5-branes wrapped around the entire compact product space. The 
near horizon geometry is AdSa x x and one can formulate a two-dimensional 
CFT at the conformal boundary of the AdSs. This CFT has central charge c = 6 NiN^ 
and A/" = (4,4) supersymmetry. (For reviews, see [46,47].) 

The moduli space of the CFT is twenty-dimensional and includes a special point 
known as the “orbifold point.” In the torus case, the theory at this point is a symmetric 
orbifold CFT with target space This is analogous to free super Yang- 

Mills theory in the AdS 5 /CFT 4 duality, and the CFT at this point is Y = Y 1 Y 5 copies 
of a c = 6 free CFT (the seed theory), composed of 4 real bosons and their fermionic 
superpartners, with an S'at orbifold permuting the N copies. Though built out of free 
theories, the symmetric orbifold structure couples the various conformal families of the 
seed theories so that at large N the quasiparticle picture of entanglement propagation 
is qualitatively altered, as we see in the following subsection. 


Free theory 

We will focus on the second Renyi entropy of the offset intervals in the thermal double, 
since, as discussed in section 3.3, this is hxed by the torus partition function. The 
spectrum of the orbifold is known exactly [48], and in appendix A, we extract from 
this the vacuum contribution to the torus partition function as g —)■ 0 , g —?• 1 : 

?) « ( 4 . 32 ) 


The limit leading to (4.32) is taken with N held hxed. This theory has a very large 
number of conserved currents; it is current dominated and Ccurrents = c. However, the 
orbifolding leads to the 1 /Y! prefactor, which means there are far fewer currents at 
large N than in a simple product of free theories. 

In the scaling limit /3 —)■ 0 with N held hxed, we can ignore the coefficient as we did 
in section 3.3, and the orbifold theory exhibits the quasiparticle dip in the 2nd Renyi 
entropy. On the other hand, the coefficient of the quasiparticle singularity in (4.32) is 
suppressed by N\ ~ Therefore to see the quasiparticle behavior, we must take 


t L D 


> log c . 


(4.33) 


At shorter times, the large coefficient suppresses the quasiparticle dip and (4.32) may 
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not be the dominant term. 

The free CFT is believed to be dual to a string theory with vanishing string tension 
(see for example [49,50]). It would be interesting to interpret (4.33) in the string theory. 
Curiously, the time scale /Slogc is the scrambling time in the gravity limit, but here 
plays a different role. 

Deformed theory 

Semiclassical gravity is far in moduli space from the free orbifold CFT. Roughly speak¬ 
ing, the coupling constant in the CFT is an exactly marginal deformation that cor¬ 
responds to the string tension, and gravity is a good approximation when the string 
tension is very large. 

Little is known about how to follow this deformation in general, though many quan¬ 
tities protected by supersymmetry can be matched in the two limits. To understand 
the behavior of the second Renyi entropy, we are interested in the spectrum of con¬ 
served currents as a function of the coupling constant. In the free limit, it is given 
by (4.32). In the gravity limit, conserved currents correspond to massless gauge helds 
in supergravity, and the only such helds are the graviton, gravitino, and gauge helds 
required by (4,4) supersymmetry. Thus all other conserved currents contributing to 
(4.32), with dimension (0, h) at zero coupling, must pick up a non-zero left-moving 
conformal weight when we move to the gravity point in moduli space. 

It seems likely that all of these currents are lifted even at leading order in the de¬ 
formation away from the free orbifold point, since they are not protected. This was 
recently conhrmed explicitly for the low-lying currents using conformal perturbation 
theory [51]. The leading deformation, which has been studied perturbatively for ex¬ 
ample in [52-56], is by an FAr-twist operator, so it exists for any N > 2. Therefore, we 
expect the CFT dehned by a small deformation of the orbifold theory with c > 12 to 
be a tractable example of a CFT with Ccurrents < c. 

5 Conclusion 

We have argued that there is a qualitative difference in the entanglement properties of 
theories with Ccurrents = c and Ccurrents < c. Although CFTs with fewer conserved cur¬ 
rents seem more generic in some sense, most of the well known CFTs in 1-1-1 dimensions 
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are rational and therefore Ccurrents = c. What is the simplest theory with Ccurrents < c? 
We do not have an answer to this question, but will suggest two candidatesd^ 

The hrst natural candidate is a theory with several scalars, and an arbitrary quartic 
potential. If the quartic potential is 0{N) invariant then the hxed points are well 
studied, but for a generic potential, the symmetry is reduced. If such a theory has a 
hxed point with reduced symmetry then it seems likely to have Ccurrents < c. 

Another class of well known examples are supersymmetric sigma models with a 
Calabi-Yau target space. A generic Calabi-Yau has no isometries, so the sigma model 
has only the conserved currents required by superconformal symmetry. A simple ex¬ 
ample of this type is the minimal version of the D1-D5 CFT discussed in section 4.4. 
It is a deformation of the orbifold theory with target (T^)^/Z 2 . This theory has central 
charge c = 12. The A/" = (4,4) superconformal algebra consists of the stress tensor, 
four /^-currents, and four supercurrents, so assuming that other currents are lifted by 
the deformation, Ccurrents < 9 < c and it follows that the entanglement entropy differs 
from the quasiparticle picture. 

Another question left open by our analysis is how to compute SAuBit) in theories 
that have Ccurrents < c < oo. We have argued that, in the dip regime, the entanglement 
entropy differs from the quasiparticle value, since the Renyi entropies disagree with 
the quasiparticle prediction. It would be very interesting to compute the vacuum 
contribution explicitly, which is the universal answer in the holographic limit. For 
theories with no extended algebra it depends only on the central charge. 

The suppression or absence of the dip in the two-intervals entanglement entropy 
after a quench might at hrst seem to clash with holographic computations for D < L 
[10-12]. In fact, for small separation between the intervals, the holographic mutual 
information after a quench has a spike (or a bump in the case of a local quench) at 
early times, which seems in qualitative agreement with the quasiparticle picture [10-12]. 
This spike though has a different origin than the dip, and can be traced in the growth of 
Saub to saturation. For D < L, Saub starts growing as 2 t for t < D/2, and continues 
~ t till saturation at time L — D/2. When combined with the growth ~ t till t = L/2 oi 
the single interval entanglement entropies Sa, Sb, this effect produces a spike between 
D/2<t<L — D/2 in the mutual information. 

In the main text, we have discussed extensively global quenches in the boundary 
state [6-8] and thermal double model [5] , but similar analysis and conclusions hold for 

thank Leonardo Rastelli for discussions about this question. 
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the local ‘joining’ quenches discussed in [10,57]. These describe the process of two senii- 
inhnite lines joined at their endpoints at an instant of time and subsequently evolving 
as a connected, inhnite system. The computation of the two-intervals entanglement 
and Renyi entropies in these systems can be reduced to that of a four-point function of 
twist operators in the UHP and thus proceeds in an analogous way to that presented 
in section 2. 

Closely related systems are the heavy local operator quenches considered in [10, 
58-60]. These are systems quenched by the insertion of a heavy local primary with 
weights h, h ~ c that produces a localized excitation. The operator quench with 
h = h = c/32 has the same energy-momentum tensor for t > 0 as the Calabrese-Cardy 
‘joining’ quench [57] and, for intervals sufficiently distant from the quench, exhibits 
the same evolution for the entanglement entropy. Therefore, the same observations 
and conclusions apply, and, in particular, our results resolve the apparent discrepancy 
between the quasiparticle picture and the holographic calculations observed in [10]. 
Moreover, in the limit of large-c, the absence of the dip in Saub after a heavy local 
operator quench can be immediately inferred from the results of [61] for the six-point 
function of two heavy and four light operators. 

The implications of our methods might also be relevant to the light operator 
quenches studied in [25,38,62-64]. Similarly, we expect our observations on light 
cone singularities to be important in the computation of the entanglement negativity 
of two disjoint intervals A,B after a global or local quench [14,65-67], which in [14,67] 
was found to be exactly 3/4 of the quasiparticle mutual information I {A, B). 
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A Vacuum character of the symmetric orbifold 


In this appendix we calculate the vacuum character of the symmetric orbifold GFT 
with target space {T‘^)^/Sn, following [48-50], and hnd the asymptotics claimed in 
section 4.4. 

The theory on has four free bosons and four free fermions, and therefore central 
charge 6. The vacuum character of this seed theory in the R sector is 


XR{(i,y) 


n=0 £eZ 


(v 


2+v‘)n 

n=l 


(1 -^gn)2(X lgn)2 

(1 - g-)4 


(A.l) 

(A.2) 


where y corresponds to the chemical potential. The coefficients c(n, tj dehned by this 
expansion can be used to construct the generating function of the NS-sector vacuum 
characters in the orbifold theories with c = QN [49,50]: 


OO OO 

Y p"xr(<;) = n 11(1 + (_i)'p5»+'/2+i/4)-o(»,0 _ ,A,3) 

N=0 n=0 l&L 

where we have set ?/ = 1. Our goal is to find as g —)• 1. 

Dehne the NS vacuum character of the seed theory (TrNsg^'i"'^/^^), 

OO 

X{(l) = EE 

n=0 fsZ 

Mil! 

?7(r)6 

= + 4g^/2 + lOg + + 55g^ + ■ ■ ■) 


(A.4) 

(A.5) 

(A.6) 
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and the alternating character (TrNs(— 


X{q) = 

n=0 £gZ 

Mil 

T]{tY 

= - 4g^/2 + lOg - 24g3/2 ^ 55^2 ^ ^ 


(A.7) 

(A.8) 

(A.9) 


Using these we can re-express the generating fnnction (A.3) as 


N=0 


exp 


E 



E 



(A. 10) 


As g —1 , 

X{q) - x{q) -4r^ , (A.ll) 

and 

X(/) ~ -(fcr)2e'"/(2fcr) ^ 

where q' = xhe hrst eqnation in (A.ll) indicates that the seed theory has 

^currents C 6. 

We are compnting the limit g —)■ 1“ with N held hxed, so we can plug these 
asymptotics into (A. 10). The leading singularity comes from k = 1: 




(A. 13) 


This is the result (4.32) used in the discussion of the second Renyi. The character 
(A. 13) is TrNsg^°“''A"‘ in the orbifold. The full partition function Z{t,t) includes a 
sum over spin structures, but the other terms are subleading in this limit. 
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